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ABSTRACT 

We argue that rich star clusters take at least several local dynamical times to form, and so are 
quasi-equilibrium structures during their assembly. Observations supporting this conclusion include 
morphologies of star-forming clumps, momentum flux of protostellar outflows from forming clusters, 
age spreads of stars in the Orion Nebula Cluster (ONC) and other clusters, and the age of a dynamical 
ejection event from the ONC. We show that these long formation timescales are consistent with the 
expected star formation rate in turbulent gas, as recently evaluated by Krumholz & McKee. Finally, 
we discuss the implications of these timescales for star formation efficiencies, the disruption of gas by 
stellar feedback, mass segregation of stars, and the longevity of turbulence in molecular clumps. 
Subject headings: stars: formation — stars: kinematics — stars: winds, outflows 



1. INTRODUCTION 

Star clusters are born from the densest gas clumps in 
giant molecular clouds and are likely to be responsible for 
the majority of stars ever formed (Lada & Lada 2003). 
The timescale over which a clump transforms into a clus- 
ter is a basic constraint for theoretical models. If forma- 
tion takes only ~ I dynamical time, then star formation 
is a result of the global collapse of the clump that internal 
sources of feedback are insufficient to impede. If it takes 
several dynamical times, then the clump gas must reach 
an approximate equilibrium, with self-gravity resisted by 
internal sources of pressure. In this case, star formation 
is a local process within the protocluster. Runaway grav- 
itational instability occurs on scales much smaller than 
the overall cluster, involving only a small fraction of the 
mass at any given time. The cluster formation timescale 
also determines how much dynamical mass segregation 
occurs in the gas-rich phase, which may account for the 
observed central concentration of massive stars in young 
clusters (Hillenbrand & Hartmann 1998). 

Elmegreen (2000) has argued that, over a wide range 
of scales, star formation occurs in ~ 1 crossing time. For 
star clusters the observational evidence cited for short 
formation times is the age spread of stars and substruc- 
ture in their spatial distributions. We examine these ar- 
guments in detail in §2 and show that the evidence for 
rich clusters in fact points to a considerably slower forma- 
tion process. Hartmann, Ballesteros-Paredes, & Bergin 
(2001) have argued for star formation in one crossing 
time for low mass, distributed star-forming regions like 
Taurus. Tassis & Mouschovias (2004) have presented 
counter arguments. However, neither of these analyses 
apply to individual rich star clusters, on which we focus. 

We define here the important timescales for a proto- 
cluster gas clump of mass M, radius R, ID velocity dis- 



persion a, density p, and column density S. The free- 
fall time is tg = (3n/32Gp) 1 ^ 2 , and the dynamical, or 
crossing, time is i^yn = R/cr. The virial parameter, 
a vir = 5a 2 R/(GM) (Bertoldi & McKee 1992), defines 
the relationship of these two; for a v ^ ~ 1, tg s» 0.5£dyn- 
Clumps are centrally concentrated (Mueller et al. 2002), 
so the free-fall and dynamical times vary by location 
within them. We therefore define tg and tdyn as mass- 
weighted averages over the region containing 90% of the 
mass, although in practice this may be difficult to deter- 
mine due to confusion on the outskirts of clusters. Fi- 
nally, we define the formation time if or m as the time over 
which 90% of the stars in a cluster form. Note, this is 
2.3 times the exponentiation time to m the accelerating 
star formation model of Palla & Stahler (2000). 

In §2 we discuss observational evidence that rich clus- 
ters take at least several dynamical times to form, draw- 
ing on various stages of the formation process. In §3 
we use a theoretical estimate of the star formation rate 
to show that long formation times are to be expected. 
Finally, we discuss the implications of this result in §4. 

2. OBSERVATIONAL EVIDENCE FOR SLOW FORMATION 

2.1. Morphologies of Gas Clumps 

Shirley et al. (2003) observed CS ( J = 5 — > 4) emission 
from about 60 dense gas clumps containing high-mass 
star formation. Their surface densities and masses are 
very similar to more revealed rich, young star clusters, 
so they are likely to be similar objects at an early stage 
of formation. Shirley et al. determined clump aspect ra- 
tios, defined as the ratio of major and minor axes of the 
20% of peak contour, which is well detected and resolved. 
The distribution of aspect ratios peaks at 1.26 ± 0.22, 
consistent with most clumps being circularly symmet- 
ric. This morphology contrasts with simulations of rapid 
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cluster formation in which one or two gas filaments tend 
to dominate (e.g. Bate, Bonnell, & Bromm 2003). Non- 
periodic simulations with stellar feedback that run for 
longer times need to be performed to see how long it 
takes to establish spherical morphologies. From the CS 
data it is not possible to tell if the clumps of gas relaxed 
before star formation started or after. Studies of earlier 
stages of star cluster formation, as traced by the Infrared 
Dark Clouds (Egan et al. 1998), can address this. 

2.2. Morphologies of Embedded Stars 

Since molecular clumps are turbulent, stars should 
form in substructures, which will dissolve as the stars 
orbit. Thus, the amount of substructure can constrain 
the formation time. The dissipation time of a substruc- 
ture depends on whether it is bound (Scally & Clarke 
2002) and the nature of its velocity dispersion (Good- 
win & Whitworth 2004). If unbound, it dissipates in a 
time ~ idyn- Assuming all stars form in substructures at 
constant rate, this implies tf olm ~ (M/M su b)td yn , where 
Msub is the mass in substructures. If a subcluster is 
bound, rather than dissolving it will sink to the cluster 
center due to dynamical friction. The sinking time from 
radius r form is t sin k = 0.68(r form /i?) 2 (A/ In A)t dyn (Bin- 
ney & Tremaine 1987), where A is the cluster to most- 
massive-subcluster mass ratio. This assumes the clus- 
ter density varies as r~ 2 , but a different exponent would 
not substantially change the result. Assuming a constant 
star formation rate, so the last subcluster formed at least 
a time <f or m/A ago, that all stars form in subclusters, and 
that on average subclusters form at rf orm = R/2, the for- 
mation time and the mass of the largest visible subcluster 
are related by i form < 0.17(2r form /i?) 2 (A 2 / In A)t dyn . To 
improve upon these approximate analytic estimates re- 
quires global numerical simulations of stars forming from 
turbulent gas (e.g. Schmeja & Klessen 2006). 

Comparing to observations, we find that, contrary to 
Elmegreen (2000), they are consistent with long forma- 
tion times. The ONC has quite smooth contours of pro- 
jected stellar surface density with no significant substruc- 
ture (Hillenbrand & Hartmann 1998), giving no limit on 
iform (see also Scally & Clarke 2002). In IC 348, there 
are 8 small subclusters, but these contain only ^10 — 20 
stars each, roughly 20% of the total 345 stars in the clus- 
ter (Lada & Lada 1995), giving tf orm ~ 5id yn even if all 
subclusters are unbound. The largest subcluster has 18 
members, so A « 19 and tf orm < 21id yn if subclusters are 
bound. Note, we have neglected the artificial enhance- 
ment of apparent substructure by patchy extinction. If 
this effect is significant, the true ages may be even larger. 
In no case do we find rich clusters with all or most of their 
mass in substructures, which would imply if orm « i dyn . 

2.3. Momentum Flux of Protostellar Outflows 

Tan & McKee (2002) estimated the star formation 
rate in 8 clusters using the observed momentum flux of 
protostellar outflows. Outflows are magnetocentrifugally 
driven from the star and inner accretion disk (e.g. Shu 
et al. 2000), expelling a fraction f w of the mass flux 
m* reaching the star. Outflow models find that the ra- 
tio f p = f w v w /vK is constant to within ~ 30% (Najita 
& Shu 1994), where v w is the outflow velocity and vk 
is the Keplerian velocity at the equatorial radius of the 
star. Thus, the total momentum flux p w — fprii^vx is 



determined by the star formation rate and the evolution 
of protostellar radii, with only a weak dependence on the 
latter. Tan & McKee (2002) find that, for clusters with 
a Salpeter IMF from 0.1 - 120 M , p w /M* ~ 87 km 
s , where M* is the total star formation rate. Loss of 
momentum in outflow-outflow interactions would lower 
this estimate, but this effect is probably small as outflow 
jets are well-collimated. Thus, observations of p w give 
a measurement of M* and hence the cluster formation 
timescale. The measurement of p w is quite uncertain, 
but the data suggest that it would take ~ 3 — 5 dynami- 
cal times to transform ~ 30% of the mass into stars (Tan 
& McKee 2002). Although this is currently one of the 
more uncertain methods of estimating tf orm , its accuracy 
should improve as models and observations of individual 
protostellar outflows are refined. 

2.4. Age Spreads of Pre- Main- Sequence Stars 

The age spread of stars in a cluster is a direct mea- 
sure of the formation time. Ages can be difficult to de- 
termine because they require fitting observed luminosi- 
ties and temperatures to pre-main sequence models that 
have quite large systematic uncertainties, particularly for 
low (sub-solar) mass stars. The models also depend on 
uncertain parameters such as the deuterium abundance, 
the accretion rate, and accretion geometry (Stahlcr 1988; 
Palla & Stahler 1999, hereafter PS99; Hartmann 2003), 
which influence the position of the "birthline", where 
stars first appear on the HR diagram. Fortunately, the 
initial stellar contraction is quite rapid, so the birthline 
position mostly affects ages < 1 Myr. Patchy extinction 
also introduces systematic errors into age estimates, as 
do unresolved binaries, although for these one may make 
approximate statistical corrections (PS99). Photometric 
variability is another potential source of error, but recent 
observations by Burningham et al. (2005) find that vari- 
ability cannot mimic age spreads of several Myr. In sum- 
mary, age estimates for intermediate mass (M ~ Mq), 
older (age > 1 Myr) stars are reasonably robust. 

In the ONC, PS99 estimate ages of 258 stars with 
masses 0.4 < tti^/Mq < 6.0 from the sample of Hil- 
lenbrand (1997). PS99 found 82 stars aged 0-1 Myr, 57 
aged 1-2 Myr, 34 aged 2-3 Myr, 17 aged 3-4 Myr, 8 aged 
4-5 Myr, 8 aged 5-6 Myr, 8 aged 6-7 Myr, and 6 aged 
7-10 Myr. Hartmann (2003) has argued that the oldest 
ages (~ 10 Myr) may be due to a problem of foreground 
contamination. We conclude that a significant fraction 
of the stars are 3 Myr old. This is a lower limit to tf orm 
since star formation is still continuing in the cluster, and 
because the sample is potentially incomplete for the old- 
est low mass stars. This result is broadly consistent with 
ONC age determinations based on Li abundances in pre- 
main-sequence stars (Palla et al. 2005), a few of which 
imply ages as large as ~ 10 Myr. Elmegreen (2000) esti- 
mates a density of nn = 1.2 x 10 5 cm -3 for the gas clump 
from which the ONC formed, giving a dynamical time of 
2.5 x 10 5 yr and an age > 12 crossing times. Note that 
Elmegreen argued star formation was rapid in the ONC, 
adopting an age spread of only 1 Myr. 

Elmegreen's dynamical time is probably somewhat low; 
a better estimate is i dyn = 0.95(a vir Gj" 1/2 (M/S 3 ) 1 / 4 . 
We use this rather than relying on a measured velocity 
dispersion because, for star clusters, this is often hard to 
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determine due to incompleteness, confusion, and varia- 
tion in the velocity dispersion with location in the clus- 
ter. The mass and surface density are easier to measure. 
For gas systems the typical virial parameter is a v i r ~ 1.3 
(McKee & Tan 2003). For stellar systems we adopt a 
King model, which implies a v i r > 2.0 if we take a to 
be the dispersion of the Maxwellian velocity distribu- 
tion. For the ONC, Hillenbrand & Hartmann (1998) find 
M = 4600 M Q , £ = 0.12 g cm" 2 , giving t dyn = 7 x 10 5 
yr and an implied cluster age of 4 dynamical times. 

In addition to the ONC, which has the best obser- 
vational data, we can make estimates of the formation 
and dynamical times for only a few other rich, forming 
star clusters. Palla & Stahler (2000) find age spreads 
of iform « 2.3 Myr for both p Ophiuchi and IC 348. For 
p Oph, the central cluster is embedded in the L1688 dark 
cloud. The cluster has a radius of w 1 pc, and the gas 
mass within this radius (which dominates the total mass) 
is roughly 1500 Mq (Lorcn 1989). The inferred dynam- 
ical time is 7.6 x 10 5 yr, implying a formation time of 3 
dynamical times - in a system that is still gas-dominated. 
The central cluster in IC 348, which contains roughly half 
the stars and for which Palla & Stahler make their age 
estimate, has a radius of 0.5 pc and contains « 200 Mq 
of stars (Lada & Lada 1995). Stars probably dominate 
the mass (Herbig 1998), so we infer a dynamical time of 
6 x 10 5 yr, giving a formation time of 4 dynamical times. 

It is more difficult to determine both formation and 
dynamical times for other systems, so, contrary to 
Elmegreen (2000), few conclusions can be drawn. For 
example, Forbes (1996) did not find evidence for an age 
spread in NGC 6531, but the analysis was insensitive 
to timescales shorter than > 3 Myr. Hodapp & Dcanc 
(1993) determined ages up to 6 Myr for stars in L1641, 
but studied only 12 objects, did not correct for binarity, 
and used relatively old pre-main sequence tracks. Palla & 
Stahler (2000) found formation times for Taurus- Auriga, 
Lupus, Chamaeleon, Upper Scorpius, and NGC 2264 
that are > 3 Myr, but the first three of these are not rich 
clusters, Upper Scorpius has undergone too much dy- 
namical spreading to reliably estimate its dynamical time 
at formation, and NGC 2264 is too distant for a complete 
census to reliably determine its mass and column density, 
and thus its dynamical time. In NGC 3603, Eisenhaucr 
et al. (1998) find both young stars (< 0.5 Myr) and 
Wolf-Rayet stars 2-3 Myr old, but as with NGC 2264 
it is too distant to reliably determine a dynamical time. 
Elmegreen (2000) comments that star clusters with age 
spreads - 10 Myr, e.g. NGC 1850, NGC 2004, NGC 
4755, NGC 6611 and the Pleiades, could be the result of 
"multiple and independent star formation events" , which 
then form a single cluster by merging or only appear to 
be a single cluster because of projection. We would argue 
that in the former case, the age spread is a true indication 
of the cluster formation time, although, as Elmegreen 
points out, the relevant dynamical timescale may need 
to be evaluated over a larger region that hosted the ini- 
tial subclustcrs. However, in all these systems, there is 
no evidence to favor merging of independent subclusters 
over continuous formation in situ. 

2.5. Age of a Dynamical Ejection Event in Orion 

Dynamical ejection events provide another method of 
age estimation. One such event involving 4 massive stars 



(a binary and two singles) that appear to have come from 
the ONC has been dated to ~ 2.5 Myr ago (Hoogerw- 
crf, dc Bruijnc & de Zeeuw 2001). The central value 
of the time since this ejection event is 2.3 ± 0.2 Myr 
in Hoogerwerf et al.'s analysis; however, if the cluster's 
distance of about 450 pc is adopted, then the best es- 
timate is 2.5 Myr. The identification with the ONC is 
based upon the extrapolation of the motion of the center 
of mass of the four stars from the ejection event to the 
present day, leading to a predicted position coincident 
with the ONC (uncertainties are a couple of pc). This 
result implies that 2.5 Myr ago the ONC was already a 
rich cluster containing at least four stars of spectral type 
earlier than 09/BO. Before this the stars had to form 
and have enough time to find and eject each other in a 
close interaction. Thus the estimate of 2.5 Myr is again 
a lower limit to tf olm for the ONC, so that tf olm > 4id yn - 

3. THEORETICAL FORMATION TIMESCALE 

Krumholz & McKee (2005, hereafter KM05) estimate 
the star formation rate in supersonically turbulent gas, 
and we use this result to compute how long star for- 
mation in a clump must continue to reach ~ 30% effi- 
ciency. Consider a clump with density and pressure pro- 
files p oc r~ kp and P oc r~ kp . Hydrostatic equilibrium 
requires that its mass, radius, and effective sound speed 
c = (P/p) 1/2 be related by M = k P c 2 R/G. The effective 
sound speed is related to the ID velocity dispersion a c \ 

by c = <i>g 2 a c \, where <frg is a factor accounting for mag- 
netic support. The clump velocity dispersion in terms of 
clump mass and surface density is thus 

-c 1 =(^P) 1/4 ^2.4M 3 1 / 4 ^ 4 kms-\ (1) 

where M 3 is the clump mass in units of 10 3 M , S 
is column density in g cm -2 , and the numerical eval- 
uation is for the fiducial parameters of McKee & Tan 
(2003) (k P = 1, 4> B = 2.8, and a vir = 1.3). Based on 
analysis of the structure of turbulent media and com- 
parison to numerical simulations, KM05 find that the 
fraction of the mass forming stars every free-fall time is 
SFR ff w 0.073 a-° T 6S M-°- 32 , where M = a cX /c s is the 
Mach number, and c s is the thermal sound speed. For 
our clump, this gives SFR ff _ cl « 0.027 T 1 - 16 (M 3 So)-° 08 , 
where T\ is temperature in units of 10 K. At this rate, 
turning 30% of the gas into stars takes 5 — 6 dynamical 
times for a 1000 Mq clump, with a very weak dependence 
on temperature, mass, or surface density. 

KM05 modeled clouds that were not centrally concen- 
trated, so we can improve this estimate by considering 
density variation. In turbulent media that are not cen- 
trally condensed, the velocity dispersion increases with 
length scale I as a oc I 1 / 2 regardless of position in the 
medium, but the tidal field of a clump may introduce a 
dependence on the distance r from the clump center as 
well. The expected variation is a oc r 1 ~ k p/ 2 : and obser- 
vations show k p ~ 1.5 (Mueller et al. 2002). On length 
scales I -C r, the tidal field of the clump is negligible and 
we should find a oc I 1 / 2 , as for a uniform medium, while 
for i 3> r the tidal field will dominate. Since the size scale 
of a star-forming core is much less than r over the vast 
majority of a star-forming clump, we can approximate 
this behavior by a « acXr 0--k p )/2 R k p /2~i t i/2 for £ ^ r 
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Thus, the star formation rate varies within the clump as 
SFR ff = SFRff-dr - 32 ^/ 2 - 1 )^ - 3 ^ 1 -^/ 2 ). Similarly, 
the free-fall time as a function of distance from the clump 
center is t$ = f ff _ cl [3/(3-fc p )] 1 / 2 (r/i? cl ) fc p/ 2 . Integrating 
over radial shells, the total star formation rate is 



471T p 



SFR ff 
ts 



dr 



(3-fc P ) 3/2 



2.3(2- k p ) 



M 



SFR, 



ff-cl 



(2) 

Thus, central condensation increases the star formation 
rate by a factor of (3 — fc p ) 3 / 2 / [2.3(2 — k p )] relative to a 
uniform medium. For our fiducial k p — 1.5, this means 
that the time required to reach 30% star formation effi- 
ciency is reduced by a factor of 1.6 relative to our previ- 
ous estimate, giving 3 — 4 dynamical times. 

4. DISCUSSION 

We have presented observational and theoretical evi- 
dence that rich star clusters require at least several dy- 
namical timescales to form. This is significant, because 
it implies that star clusters cannot form by a process 
of freely decaying turbulence leading to global collapse, 
which could not possibly take so long. Instead, some- 
thing must impede or entirely prevent global collapse, 
so that rich star clusters are in approximate equilibrium 
during their assembly. For the ONC, p Oph, and IC 348, 
where we can reasonably estimate both formation and 
dynamical times, formation typically requires > 3 — 4 
dynamical times, consistent with theoretical predictions 
for the time required to turn ~ 30% of the gas into stars. 
The central regions have much smaller dynamical times 
than the cluster average, so we predict that in the center 
a larger fraction of the gas will form stars. 

Another implication of this work is that massive star 
feedback may not be as effective as once assumed in dis- 
persing gas in young clusters. This is consistent with 
observations that massive stars are not always the last 
to form in their clusters (e.g. Eisenhauer et al. 1998; 
Hoogerwerf ct al. 2001), and theoretical work showing 
that dumpiness greatly inhibits gas dispersal (Tan & Mc- 
Kee 2004; Dale et al. 2005). Tan & McKee find that, in 
a clump like the proto-ONC, if ~ 3% of the mass forms 
stars per dynamical time, feedback requires ~ 2 Myr 
(about 3 tdyn) to disperse the gas. Dynamical ejection 



of massive stars, as observed in the ONC (Hoogerwerf et 
al. 2001; Tan 2004), would increase this time. 

Long formation times are also important for mass seg- 
regation. For example, a 3 Myr formation time for the 
ONC corresponds to about 8 diameter crossing times at 
the half mass radius of 0.5 pc, which is the unit of time 
used in the study of Bonnell & Davies (1998). If the 30 
most massive stars are born at the half-mass radius then 
after 3 Myr the median location of the 6 most massive 
stars migrates to only 0.075 pc, suggesting that some 
of the observed mass segregation (Hillenbrand & Hart- 
mann 1998) could be dynamical rather than primordial. 
Gas drag will likely enhance the segregation beyond the 
purely N-body effects explored by Bonnell & Davies. 

Although undriven supersonic turbulence decays in 
~ 1 dynamical time (Stone, Ostriker, & Gammie 1998), 
the observed turbulence in molecular clumps does not 
damp so quickly. Driving by protostellar feedback is 
a possible explanation. A virialized clump that radi- 
ates away half its kinetic energy per dynamical time 
has a luminosity L diss w 4M 3 5/4 Eo' 4 L Q , but the IMF- 
averaged energy release associated with accretion is 

L acc » 2000M 3 £o /4 (77/TO)- 1 L Q , where T] is the num- 
ber of dynamical times required to transform 50% of the 
gas into stars. Outflows should eject about half this en- 
ergy back into the clump (Shu et al. 2000), so even if 
only ~ 1% of this goes into driving turbulence, that is 
sufficient to offset the decay. Recent observations that 
find outflows inject enough energy to maintain turbu- 
lence (Williams, Plambeck, & Heyer 2003; Quillcn ct al. 
2005) support this idea, as do the numerical simulations 
of Li & Nakamura (2006). However, this work is prelim- 
inary and has not yet shown that feedback can maintain 
turbulence over a cluster lifetime of 4 td yn that observa- 
tions seem to require. 
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